Spin-polarized electron transport in ferromagnet/semiconductor heterostructures: 
Unification of ballistic and diffusive transport 
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A theory of spin-polarized electron transport in ferromagnet/semiconductor heterostructures, 
based on a unified semiclassical description of ballistic and diffusive transport in semiconductor 
structures, is developed. The aim is to provide a comprehensive framework for systematically 
studying the interplay of spin relaxation and transport mechanism in spintronic devices. A key 
element of the unified description of transport inside a (nondegenerate) semiconductor is the ther- 
moballistic current consisting of electrons which move ballistically in the electric field arising from 
internal and external electrostatic potentials, and which are thermalized at randomly distributed 
equilibration points. The ballistic component in the unified description gives rise to discontinuities 
in the chemical potential at the boundaries of the semiconductor, which are related to the Sharvin 
interface conductance. By allowing spin relaxation to occur during the ballistic motion between the 
equilibration points, a thermoballistic spin-polarized current and density are constructed in terms 
of a spin transport function. An integral equation for this function is derived for arbitrary profile 
of the electrostatic potential and for arbitrary values of the momentum and spin relaxation lengths. 
Detailed consideration is given to field-driven spin-polarized transport in a homogeneous semicon- 
ductor. An approximation is introduced which allows one to convert the integral equation into a 
second-order differential equation that generalizes the standard spin drift-diffusion equation. The 
spin polarization in ferromagnet/semiconductor heterostructures is obtained by invoking continuity 
of the current spin polarization and matching the spin-resolved chemical potentials on the ferromag- 
net sides of the interfaces. Allowance is made for spin-selective interface resistances. Examples are 
considered which illustrate the effects of transport mechanism and electric field. 

PACS numbers: 72.25.Dc, 72.25.Hg, 73.40.Cg, 73.40.Sx 



I. INTRODUCTION 

Considerable attention has been devoted in the past 
few years to the study of spin-polarized electron trans- 
port in hybrid nanostructures composed of different types 
of material, such as nonmagnetic or magnetic semicon- 
ductors, normal metals, ferromagnets, and superconduc- 
tors. The motivation behind these efforts derives from 
the desire to understand the principles of operation, 
to assess the performance, and to explore the field of 
possible applications, of solid-state devices relying on 
the control and manipulation of electron spin ("spin- 
tronic devices" ) i2*2*4 Particular emphasis in spintron- 
ics research is currently placed on the study of spin- 
polarized transport in heterostructures formed of a non- 
magnetic semiconductor and two (metallic or semicon- 
ducting) ferromagnetic contacts ^'-' 6 ' '- Structures of this 
kind are considered promising candidates for future tech- 
nological applications. For the actual design of spin- 
tronic devices, a detailed theoretical understanding of 
spin-polarized transport in ferromagnet/semiconductor 
heterostructures is indispensable. Up to now, a number 
of pertinent studies have been undertaken, which mostly 
rely on the drift-diffusion model. 

Schmidt et al& describe the spin polarization by the 
same diffusion equation for the chemical potential as 
used to treat spin-polarized transport in ferromagnet /- 
normal-metal heteroj unctions AiMiiiS For a metallic- 
ferromagnet/semiconductor heterojunction with a per- 



fect interface (no interface resistance or spin scattering) , 
they find that, as a consequence of the large conductance 
mismatch, the injected current spin polarization is very 
low. Filip et idMi and Rashba 14 suggest that efficient spin 
injection can be obtained by introducing spin-selective 
interface resistances, for example, in the form of tunnel 
barriers. This idea is pursued in a number of detailed the- 
oretical investigations in which the interface resistances 
are taken into account either phenomenologically by in- 
troducing discontinuities in the chemical potentials at the 
interfacer^iSiiLi2iiS, or explicitly by treating the Schot- 
tky barrier arising from band bending in the interface 
depletion regioni2Si2i Yu and Flatt o 18 ' 19 have derived a 
drift-diffusion equation for the density spin polarization, 
which allows, in particular, the effect of applied electric 
fields to be studied. A formalism taking into account 
the effect of the electron-electron interaction on spin- 
polarized transport in metals and doped (degenerate) 
semiconductors in the diffusive regime has been intro- 
duced by D'Amico and Vignale, 2 • and has subsequently 2 ^ 
been generalized to include the effect of applied elec- 
tric fields. Spin injection under conditions where, in 
the semiconductor, ballistic transport prevails over drift- 
diffusion has been studied by Kravchenko and Rashba 24 
within a Boltzmann equation approach; they find that 
in the absence of spin-selective interface resistances, the 
Sharvin interface conductance^ controls the injection ef- 
ficiency. Spin injection across a Schottky barrier, arising 
from thermionic emission as well as tunneling injection, 
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has been treated by Shen at aliSi within a Monte Carlo 
model. Phase-coherent tra nsport in the ba llistic limit has 
been studied, e.g., in fiefs. I27l28l2c|3ri3ll32ll33 

It emerges that the theory of spin-polarized electron 
transport in ferromagnet/semiconductor heterostruc- 
tures has reached a level of considerable sophistication. 
Nevertheless, we believe that certain aspects of the semi- 
conductor part of the transport problem require a more 
detailed, unified treatment, such as the interplay of spin 
relaxation and transport mechanism all the way from the 
diffusive to the ballistic regime, and the effects of the spa- 
tial variation of the electrostatic potential profile. In the 
present work, we provide a comprehensive framework for 
systematically dealing with these aspects. 

The starting point is our unified semiclassical descrip- 
tion of (spinless) ballistic and diffusive electron transport 
in parallel-plane semiconductor structures^ in which the 
idea of a thermoballistic transport mechanism was intro- 
duced. The latter relies on the concept of a "thermobal- 
listic current" inside the semiconducting sample. This 
current consists of electrons which move ballistically in 
the electric field arising from internal (built-in) and ex- 
ternal electrostatic potentials, and which are thermalized 
at randomly distributed equilibration points (with mean 
distance equal to the mean free path, or momentum re- 
laxation length) due to coupling to the background of 
impurity atoms and phonons. The current-voltage char- 
acteristic for nondegenerate systems as well as the zero- 
bias conductance for degenerate systems are expressed 
in terms of a reduced resistance; for arbitrary momen- 
tum relaxation length and arbitrary potential profile, the 
latter quantity is determined from a resistance function 
which is obtained as the solution of an integral equa- 
tion. The thermoballistic chemical potential and current 
are derived from this solution as well. The chemical po- 
tential exhibits discontinuities at the boundaries of the 
semiconductor, which are related to the Sharvin interface 
conductance. 

In order to develop, within the unified description, a 
theory of spin-polarized electron transport in (nondegen- 
erate) semiconductors, we introduce the "thermoballistic 
spin-polarized current" which generalizes the thermobal- 
listic current of Ref. 34 so as to allow spin relaxation to 
take place during the ballistic motion between the equi- 
libration points. The thermoballistic spin-polarized cur- 
rent is constructed in terms of a "spin transport function" 
that determines the spin polarization inside the semicon- 
ductor for arbitrary potential profile and arbitrary val- 
ues of the momentum and spin relaxation lengths. This 
function satisfies an integral equation which follows from 
the balance equation connecting the thermoballistic spin- 
polarized current and density. The appearance of an 
integral equation in the unified description of electron 
transport (with or without account of the spin degree 
of freedom) reflects the nonlocal character of the trans- 
port across the ballistic intervals between the equilibra- 
tion points. For electron transport in a homogeneous 
semiconductor without space charge, driven by an ex- 



ternal electric field, the integral equation for the spin 
transport function can be converted, in an approxima- 
tion that is adequate for demonstrating the principal ef- 
fects of the transport mechanism, into a second-order 
differential equation that generalizes the standard spin 
drift-diffusion equation. The spin polarization along a 
ferromagnet/semiconductor heterostructure is obtained 
by invoking continuity of the current spin polarization 
at the interfaces and matching the spin-resolved chemi- 
cal potentials on the ferromagnet sides of the latter, with 
allowance for spin-selective interface resistances. 

As a prerequisite to developing a theory of spin- 
polarized electron transport in semiconductors within the 
unified transport model of Ref. I34I we have to modify 
and complete the formulation given in that reference. 
This will be done in the next section. In Sec. Ill, the 
spin degree of freedom is introduced into the unified de- 
scription. The integral equation for the spin transport 
function inside a semiconducting sample is derived, and 
the generalized spin drift-diffusion equation is obtained. 
Spin-polarized transport in heterostructures formed of a 
nonmagnetic, homogeneous semiconductor and two fer- 
romagnetic contacts is treated in Sec. IV. We demon- 
strate the procedure for the calculation of the current 
and density spin polarizations across a heterostructure in 
the zero-bias limit and of the injected spin polarizations 
for field-driven transport. Various examples are consid- 
ered which illustrate the effects of transport mechanism 
and electric field and exhibit the relation of the unified 
description to previous descriptions by other authors. In 
Sec. V, the contents of the paper are summarized and 
our conclusions as well as an outlook towards applica- 
tions and extensions of the present work are presented. 
In the Appendix, details of the extended formulation of 
the unified transport model outlined in Sec. II are worked 
out. 



II. UNIFICATION OF BALLISTIC AND 
DIFFUSIVE TRANSPORT IN 
SEMICONDUCTORS 

The unified description of (spinless) ballistic and diffu- 
sive electron transport developed in Ref. [34| has yielded, 
for the nondegenerate case, the current-voltage charac- 
teristic for a semiconducting sample enclosed between 
two plane-parallel contacts. There, the discontinuity of 
the chemical potential has been placed at the interface 
at one or the other end of the sample; this gave rise to 
an ambiguity in the behavior of the chemical potential 
inside the sample. In order, nevertheless, to obtain a 
unique current-voltage characteristic, the reduced resis- 
tance determining the latter was subjected to a heuristic 
symmetrization procedure (see Sec. IV. C of Ref. 0). 

In the following, we extend in a systematic way the 
unified description by treating simultaneously the ef- 
fects of the two interfaces on an equal footing. In 
this way, we arrive at unique chemical potentials, cur- 
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FIG. 1: Schematic diagram showing a semiconducting sam- 
ple of length S enclosed between two plane-parallel ferro- 
magnetic contacts. Illustrated are expression 12. 161 for the 
thermoballistic current J(x) and expression 13.281 for the off- 
equilibrium thermoballistic spin-polarized current J-(x). 



rents, and densities inside as well as at the ends of 
the semiconducting sample. This is prerequisite to the 
study of spin-polarized electron transport in ferromag- 
net/semicinductor heterostructures, which is the princi- 
pal aim of the present work. As in Ref. 0, we work 
within the semiclassical approximation, ignoring all co- 
herence effects. 



A. Thermoballistic transport 

We consider a semiconducting sample bordering on a 
left contact at x = x\ and on a right contact at x = X2 
(see Fig. 0, so that S = x^ — x\ is the sample length. 
The geometry of the set-up is one-dimensional, whereas 
the electron motion is treated in three-dimensional space. 
Reformulating the unified description of electron trans- 
port in a nondegenerate semiconductor^ we introduce 
the (net) electron current density (electron current, for 
short) J(x',x") across the ballistic interval [x',x"] be- 
tween two equilibration points x' and x" , 



J(x',x") = v e N c e-P E ™( x '< x ") 



e /3n(x') _ e f3fj.(x") 



(2-1) 



(xi < x' < x" < X2), which is the difference of the 
ballistic current J l (x',x") injected at the left end at x' , 



(2.2) 



and the analogous ballistic current J r (x',x") injected at 
the right end at x" , 



J r (x',x") = v e N c e-^ E ^ x '' x "^ x "^ 



(2.3) 



Here, the function n(x) is the chemical potential at 
the equilibration point x. [In Ref. l34l the term quasi- 



Fermi energy (symbol Ep) was used for the chemi- 
cal potential as defined, e.g., in Ref. |35[] Further- 
more, v e = (27TTO*/?)- 1 / 2 is the emission velocity, N c — 
2(27rm*//9/i 2 ) 3 / 2 is the effective density of states at the 
conduction band edge, to* is the effective mass of the elec- 
trons, and (i = {ksT)- 1 . The currents (£3 and 
contain only the electrons transmitted across the sample, 
i.e., the electrons with sufficient energy to surmount the 
potential barrier determined by 



E™(x',x") = E™(x',x")-E°; 



(2.4) 



here, E™(x' ,x") is the maximum value of the potential 
profile E c (x) in the interval [x',x"], and E® is its overall 
minimum across the sample. The profile E c {x) is, in 
general, a self-consistent solution of a nonlinear Poisson 
equation involving the conduction band edge potential 
and the external electrostatic potential. 

Expressions (|2.2|) and l|2.3|l for the currents injected at 
the left and right ends of a ballistic interval follow from 
Eqs. (20) and (21) of Ref. [3 if classical transmission 
probabilities are used. Tunneling effects can be included 
by replacing the classical probabilities with their quan- 
tal (WKB) analogues, as done, e.g., in Ref. |3^. In the 
present paper, we do not consider this possibility. 

It is convenient to rewrite Eq. I|2.1|l in the form 



J(x',x") 
the quantity 



[J{x')-J{x")] 



J(x') = v e N c e~^ E °-^ 



(2.5) 



(2.6) 



is the current injected at the left end x' of the ballistic 
interval into the right direction in a flat profile E c (x) = 
E®, and similarly for the current Jix") injected at the 
right end x" into the left direction. 

The ballistic density, at position x in the interval 
[x',x"], of transmitted electrons injected at the left end 
at x' is given by 

n l (x',x";x) = ^-{2(3/ vm*) 1 1 2 

x / dpe~ p[p2/2m * +Eci - x) ~ li(x ' )] 
Jo 

x 6(p 2 /2to* + E c (x) - E™(x', x")) 
= ^lC(x',x";x)e-^^'' x "1-^ x '^ , 

(2.7) 

where 

C(x',x";x) = e ^(*'^")~E c { x )\ 

x C AcmE?{x' ,x") - E^x)]) 1 ' 2 ) ■ 

(2.8) 

the function erfc(a;) is the complementary error func- 
tion^ and < C(x\ x"; x) < 1 [the conserved currents 
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(|2.2() and l|2.3[) are, of course, independent of x]. Anal- 
ogously, the ballistic electron density in the transmitted 
current injected at the right end at x" is 

n r (x',x";x) = ^C(x',x";x)e-^ E ™( X '< X "^^ X "V . 

(2.9) 

Dividing the current by the density [or 

by Q2.9|l ]. we obtain for the average velocity v(x',x";x) 
at position x of the electrons injected at either end of, 
and transmitted across, the interval [x 1 ,x"] 



v(x', x"\ x) 



2v R 



J l ' r {x',x") _ 

n l ' r (x',x";x) C(x',x";x) 



(2.10) 



which depends only on the potential profile. For constant 
profile, one has C(x',x";x) = 1, and the electrons move 
with the average velocity 2v e from one or the other end to 
its opposite. This is the average velocity of the injected 
electrons also in the case of a position-dependent profile. 
However, some of these electrons are reflected back, so 
that the average velocity at position x of those electrons 
which have sufficient energy to pass over the top of the 
potential profile must be higher than 2v e , namely, equal 
to the velocity given by Eq. (|2.1U|) . 

In analogy to Eq. 112. ip . the sum of the densi- 
ties n(x',x";x) and n r (x', x"; x) is the ballistic den- 
sity n(x' ,x" ;x) of transmitted electrons in the interval 
[x',x"}, 



n(x' , x"; x) 



C{x' 1 x"]x)e-~ f3E ' {x '- x "^ 



(2.11) 



The density n(x', x"; x), like the current J(x',x"), com- 
prises only electrons that participate in the transport. 

From the ballistic current (|2.5|) . the thermoballistic 
current J(x) at position x inside the semiconductor is 
constructed by summing up the weighted contributions 
of the ballistic intervals \x\ x"\ for which x\ < x' < x < 
x" < x 2 [see Eq. (23) of Ref.H|, 

J(x) = w-l(xx,x 2 ;1) [Ji - J 2 ] 

r x dx' 



1 



■w 2 (x',x 2 ;l)[J(x') - J 2 ] 



Xl 

X2 dx' 



w 2 {x u x';l) [Ji-J(x')] 



f x dr' f X2 dx" 

(2.12) 



with 



w n (x',x";l)=p n (\x"-x'\/l)e-P E ?( x '< x ") (2.13) 

(n = 0,1,2,3), where I is the momentum relaxation 
length of the electrons, which comprises the effect of re- 
laxation due to electron scattering by impurity atoms 



and phonons. The probabilities p n (x/l) of occurrence 
of the ballistic intervals depend on the dimensionality of 
the transport [note that w n (x',x";l) is symmetric with 
respect to an interchange of x' and x"\. In Eq. H2.12fl . 
the quantities J\ t2 = J(x\. 2 ) are fixed by the chemical 
potentials /ix,2 — ^{ x i.2) on the contact sides of the in- 
terfaces at X\ 2, i.e., immediately outside of the sample, 



1.2 



= v P N r e~^ E °-^) 



(2.14) 



For later convenience, we introduce a symbolic operator 
W(x' ,x";l) to write expression 12.12|) in the condensed 
form 

J{x) = / T" P T" w(x '' x " ] l) [J{x ' ] ~ J{x " )] ' 

(2.15) 

which, in view of Eqs. (|2.ti|) and (|2.14|l . may also be writ- 
ten as 



J(x) 



v P N c e~^ 



d^_ 

1 1 I 



X2 dx" 



(x',x";l) 

(2.16) 



Analogously, we introduce the thermoballistic density in- 
side the semiconductor, n(x), as 



n(x) 



where 



Nc e -,3E° 



x dx' f X2 dx" „ , . .. , x 
-r- / —W c (x',x";l;x) 

xi ' Jx J 



(2.17) 



W c (x', x";l; x) = W(x' , x"; I) C(x', x";x) . (2.18) 

Again, the current J(x) and the density n{x) comprise 
only electrons that participate in the transport. 

The thermoballistic current Ij2.12|l by itself is not, in 
general, conserved, but together with the background 
current J back (x) it adds up to the conserved physical 
current J: 



J(x) + J back {x) = J = const. 



(2.19) 



The background current is confined within the sample 
and, therefore, must vanish when integrated over the lat- 
ter, which implies that the thermoballistic current J(x) 
averaged over the sample yields the physical current J, 



■f'2 



1 f X2 

/ dx J(x) = J 

- xi J X1 



(2.20) 



The non-conservation of the thermoballistic current can 
be viewed as arising from a source term Q(x) associated 
with the gain of thermoballistic electron density due to 
the coupling to the background, as expressed by the equa- 
tion 



dJ(x) 
dx 



Q(x) . 



(2.21) 
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In the background, the source term appears as a sink 
term describing the loss of electron density, 



-\back 



(x) = -Q(x) 



(2.22) 



Again, since the background electrons are confined to the 
sample, the integral of Q back (x) over the sample must 
vanish and, therefore, also that of the thermoballistic 
source term Q(x), 

px 2 

/ dxQ(x)=0. (2.23) 
Owing to Eq. I|2.21[) . this implies 



J(xi) = J{x 2 ) = kJ 



(2.24) 



that is, the thermoballistic current entering at one end 
of the sample, x — must be the same as the one 



leaving at the other end, x 



The quantity k has 



been introduced for later convenience; it normalizes the 
thermoballistic current on the sample sides of the fer- 
romagnet/semiconductor interfaces to the physical cur- 
rent J. We remark that, in analogy to the thermoballis- 
tic current, the thermoballistic density as well as other 
thermoballistic quantities introduced later in the devel- 
opment all have their background complement. 

Condition (|2.2Q(I has been used in Ref. to obtain 
the fundamental integral equation for the determination 
of the thermoballistic current. Condition (|2.24() is new, 
and provides us with an extension of the formalism of 
Ref. |34| which allows us to establish a unique thermobal- 
listic chemical potential inside the sample, as will be 
shown in the following. 



B. Thermoballistic chemical potential, current, and 
density 

Substituting expression l|2.12[) in condition (|2.20() . we 

obtain 



x 2 - Xl 



J = 



X2 dx' 

ui(x ll x 2 ;l)+ I —u 2 (xi,x';l) 
X2 dx 1 

u 1 (x 1 ,x 2 ;l)+ I —ru 2 (x',x 2 - 1 l) 



X2 dd_ 
~T 



Xl 



I 



u 2 (x',x 2 ; I) - u 2 (xi,x'; I) 



Ji 



f X2 dx" 

/ ^-u 3 (x',x";l) 

J x x ' 



J(x') , (2.25) 



where 



■ w n (x', x"; I) 



(2.26) 



[note that u n (x',x";l) is antisymmetric with respect to 
an interchange of x' and x"]. Equation l|2.25|l is a basic 
condition on the function J{x), and hence, via Eq. I|2.6|l . 



on the chemical potential fJ.(x), whose determination al- 
lows all relevant transport quantities to be obtained. 

For given current J, only the value of the chemical 
potential at one of the interfaces with the contacts can 
be prescribed. If, at the interface at x\, we prescribe the 
value of /xi, i.e., that of J\ [case (i)], we can hnd the value 
of J(x) at the other interface at x 2 by re-expressing in 
Eq. (|2.25|l J 2 as J{x 2 ) and replacing x 2 with the variable 
x, thereby obtaining an integral equation for the function 
Jix) in the range X\ < x < x 2 . If, now, J{x 2 ) is required 
to assume a ^reassigned value for which we re-introduce 
the symbol J 2 , then the current J on the left-hand side 
of Eq. (|2.25() is fixed at some value J\ . We denote the 
associated solution of the integral equation by J\{x). 

On the other hand, prescribing the value /j, 2 for the 
chemical potential at the interface at x 2 [case (ii)], we 
re-express in Eq. (|2.25(l J\ as J{x\). Then, replacing X\ 
with the variable x, we obtain an integral equation for 
the function J(x) in the range x\ < x < x 2 . With Jix\) 
required to assume a preassigned value J\ , the current J 
is now fixed at the value J 2 . The associated solution of 
the integral equation is denoted by J 2 {x). 

To determine the functions J\, 2 ix) explicitly, we pro- 
ceed in the following way. In case (i), we define the "re- 
sistance function"— 

Xi(x) = Jl ~f l(x) ; Xi(zi) = 0, (2.27) 

and obtain from Eq. l|2.25|l . following the procedure out- 
lined above, 



x — X\ 



f x dx' 

ui(x\,x;l) + / — u 2 {x',x; I) 
x dx[ 

I 



■ii 



u 2 (x',x;l) - u 2 {xi,x']l) 

xi(O = 0, 



r dx" 

/ ^- U3 (x',x";l) 

J xi ' 



(2.28) 



which is a linear, Volterra-type integral equation for 
Xi{x)- Letting x — + xf in Eq. 12.28fl . we find, using 
the properties of u n {x' ', x"; I), 



xi(4) 



Ji 



(2.29) 



With this discontinuity incorporated in it, the solution 
X\{x) is unique and continuous for x\ < x < x 2 . The 
solution of Eq. lf2~2l^|l can be obtained in closed form un- 
der special conditions; in the general case, one solves 
this equation efficiently by discretization and numeri- 
cal propagation, using the initial value Xi(^i") given by 
Eq. J22<| . 

In case (ii), we define the resistance function 



X2(x) 



J2 



X2(x 2 ) = , 



(2.30) 
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which satisfies the integral equation 

x 2 — x 



I 



Ui(x, x 2 \l) + I u 2 (x,x';l) 



I 



X2(x) 



X2 <M_ 

X2 dx" 

~1 



u 2 (x',x 2 ;l) - u 2 (x,x';l) 



u 3 (x ,x ;l) 

(2.31) 

The solution Xi{ x ) is discontinuous at x = x 2 : 

X2{x -) = ^) -J-2) = ef3[ E c{x2) -E° c] X2{X2) 

(2.32) 

It follows from Eqs. I|2.28|l and (|2.31|) . using the proper- 
ties of u n (x',x"; I), that the functions Xi( x ) an d X^ix) 
are related by 



X2{x) = X*i( x o - x) 



(2.33) 



where xq — xi + x 2 ; the asterisk attached to xi indicates 
that this function is to be calculated using the reverse 
of the profile E c (x), given by E*(x) = E c (x Q — x). If 
the profile is symmetric, E*{x) = E c (x), the functions 
Xi{x) and X^{x) are the reverse of one another, X2{x) = 
Xi{x - x). 

The two functions J7i_ 2 (x) are not, in general, equal 
and yield different chemical potentials ^1.2(2;) via Eq. 

Then, in view of Eq. (|^TI|) . Eq. QTZfy implies 
/ii(xi) 7^ fix, and the chemical potential fJ.i(x) is dis- 
continuous at the interface at x = Xi, i.e., its value on 
the semiconductor side of the interface is not equal to its 
value at the interface itself. Analogously, fi 2 (x 2 ) 7^ [i 2 . 
The ambiguity thus found is a generalization of the am- 
biguity of the chemical potential in the ballistic limit 
l/S — > 00^ where it may either be associated with the 
current injected at x 2 , in which case it is discontinuous 
at x = x\ (Sharvin resistance at the interface at x = x\), 
or with the current injected at x = x\, in which case 
it is discontinuous at x = x 2 (Sharvin resistance at the 
interface at x = x 2 ). In the Appendix, details of the con- 
struction of a unique thermoballistic chemical potential 
fj-(x), current J(x), and density n(x) in terms of the two 
solutions Xi,2(x) are presented. Here, we only summarize 
the results. 

A quantity x is introduced as 

X = aixi(x 2 ) + 02X2(2:1) , (2.34) 

where di + a 2 = 1- The coefficients a\ and a 2 are given 
by 

«1,2 

j 

a 



a\,2 



(2.35) 



where 



a 1 



dx' 



{w 2 {x U x']l) [X2{x') ~ X2{xi)] 



+ w 2 (x',x 2 ; l)x2(x')} 



(2.36) 



a 2 



X2 dx' 

— {w 2 (x 1 ,x';l)xi(x') 
11 ' 

+ w 2 (x',x 2 ;l)[ Xl (x')- X i(x 2 )]} , (2.37) 



and 



a = a\ + a 2 ; 



(2.38) 



for a symmetric potential profile E c (x), we have d\ 
d 2 = 1/2. 

The current J{x) is expressed as 



J(x) = l(Ji + J 2 )-JX-(x), (2.39) 



where 
X-(x) 



a i 



Xi(x) - 2X1(^2) 



«2 



X2{x) - ~X2{xi) 



(2.40) 



(xi < x < x 2 ). The currents Cf\ and J 2 satisfy the 
relation 



Jl-J 2 = Jx 



(2.41) 



With the use of Eq. 12.14fl . the current- voltage charac- 
teristic is then obtained in the form 



J = v e N c e~ pE ' - (1 
X 



where 



V = 



Mi - M2 



(2.42) 



(2.43) 



is the voltage bias between the contacts, and E p = 
E™{x\, x 2 ) — fx±; the quantity x, given by 



X 



e -!3E™(x u x 2 ) 



(2.44) 



is the "reduced resistance" i^i It replaces, in the present 
extended unified description, expression (58) of Ref. 0, 
which was obtained, in a heuristic way, by taking the 
mean value of the reduced resistances corresponding to 
case (i) and (ii), respectively. To determine x m the 
diffusive and ballistic regimes, we evaluate the functions 
Xi,2(a;) by following the development given in the Ap- 
pendix of Ref. 34. In the diffusive regime l/S <C 1, we 
find Xi( x 2) = X2(a^i) = X; which leads to 



1 



X 



2po(0) J X1 I 



C2 e -^[E^(xi,x 2 



-Ec(x)] 



(2.45) 



where the values of po(0) for one-, two-, and three- 
dimensional transport are given in Sec. II of Ref. In 
the ballistic limit l/S — + 00, we have x = 1. 

According to Eqs. (E3, (|zH3|) . (|!0§|) . and the 

thermoballistic chemical potential fJ.(x) is given by 





' X-{x)~ 




'1 X-{x)~ 




[2 X . 




[2 X J 



X-(x) 



(2.46) 
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(xi < x < x 2 ), where 



where 



^=2 16 



/3M1 -I- e&P- 



ie^ 2 ) 



(2.47) 



and the thermoballistic equilibrium electron density n(x) 

by 



l{x) = N ce -P[Ec(x)-»(.x)] 



(2.48) 



The thermoballistic chemical potential /x(x) [and hence 
the thermoballistic equilibrium density n{x)\ are discon- 
tinuous at the interfaces at x\ 2 , 



/?e 2 J 



1 = a 2 



Si 

/3e 2 J 



(2.49) 



(2.50) 



as can be shown with the help of Eqs. (|2.29() and l|2.32[l . 
respectively, and Eqs. (|2.4U|) and I2.42[l . Here, 

Gi,2 = I3e 2 v e n h2 (2.51) 

are the respective Sharvin interface conductances) 2 ^^ 
with 711,2 = n(xx t 2). The discontinuities of the functions 
exp[(3/i(x)] and h(x) are proportional to J. In the dif- 
fusive regime l/S <C 1, where, according to Eqs. (|2.42|) 
and (|2.45|) . J oc l/S, the discontinuities approach zero 
together with l/S. 

The thermoballistic current J(x) is obtained in terms 
of X- (x) and x by substituting expression l|2.39[l in com- 
bination with Eq. (|2.41() in Eq. I|2.15[l [or, more explicitly, 
in Eq. (E33)], 

J(x) = J \ w 1 (x 1 ,x 2 ;l)x 



+ j -^-w 2 (x',x 2 ;l) ~-X-(x') 
^-w 2 (xi,x';l) ~+X-(x') 

+ jf ^ £ ^fw 3 (x',x";l)[x-(x")- X -^)]} ■ 

(2.52) 

The thermoballistic density n(x) is found in a similar 
fashion from Eq. I|2.17|l . 

n(x) = — \x coihifieV /2)W(x 1 ,x 2 \l]x) 
2v e { 

-^-to 2 (x',x 2 ;l;x) ^+x-(x') 



XI 

X2 dx' 



— to 2 {xux\l;x) y- -x-(x) 



x dx' rX2 ^" 



dx' 1 



W 3 (x',x";i,x) [x-(x") +x-{x')]j , 

(2.53) 



W(xi,x 2 ; I; x) = XOx(x\,x 2 ; l;x) 

rX dx' 

— tt> 2 (x',x 2 ;l;x) 

rX2 dx ' < > 1 ^ 



and 



r*te r^L nai3 j^. l;x) (2 . 54) 

Jxx ' Jx ' 



to n (x', x"; I; x) = w n (x', x"; I) C(x', x"; x) . (2.55) 



In the zero-bias limit V — > 0, expression l|2.53[l reduces 
to the form n(x) = fix W(x±, x 2 ; I; x), from which the 
physical meaning of the function 2U(iei, x 2 ;l;x) becomes 
apparent. 

In the diffusive regime l/S <C 1, we have 



e W) = v+ _ Ic{xi,x) -I c {x,x 2 ) ^ 



I c (xi,X 2 ) 



where 



I C {X',X") 



(2.56) 



(2.57) 



and 



J{x) = J 



v _ du(x) 
— n(x) — — 
e dx 



dE c (x) 1 dn(x) 
n{x) — ; 1- — 



dx 



13 dx 



(2.58) 

where v = 2ev e (3{l) is the electron mobility, and (I) = 
Po{0)l is the effective momentum relaxation length. 34 
Integrating Eq. (|2.58|) . we retrieve the current-voltage 
characteristic l|2.42|l with x given by Eq. (|2.45|l . The 
thermoballistic density n(x) becomes equal to the equi- 
librium density n(x) given by Eq. (|2.48() . 
In the ballistic limit l/S — > 00, we have 



e /Jf*("0 = v+ 



2 e -PE™(xi,x 3 ) 



- (ai - a 2 )}?7- , 
(2.59) 



where d\ and a 2 are to be calculated from Eqs. (|2.35(1 - 
(j2~3%|) with 



and 



/3E™(x u x) 



{3E™(x,x 2 ) 



(2.60) 



(2.61) 
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Further, 

J(x) = J = 2v e N c e -0z™{xux2) n _ (2.62) 

and 

n(x) = N c C(x 1 ,x 2 ;x)e- (3E ™ (xi < X2) r] + . (2.63) 

As our model demands, expressions (|2.62|l and 12.63JI are 
identical to the original ballistic expressions (|2.5|) and 
(|2.11() with x 1 = x\ and x" = x 2 - 

C. Field-driven transport in a homogeneous 
semiconductor 

To illustrate the formalism developed so far, we now 
consider electron transport in a homogeneous semicon- 
ductor without space charge. The electrons are assumed 
to be driven by a (constant) electric held £ directed an- 
tiparallel to the x-axis, i.e., they are moving in a linearly 
falling potential of the form 

E c (x) = E c { Xl ) - e\£\(x - Xl ) , (2.64) 

in which case 

C(x',x";x) = C(e(x-x')) 

= e eix ~ x,) erfc([e{x - x')} 1/2 ) , 

(2.65) 

where e = /?e|£|. The latter quantity is related to the 
voltage bias V via eS — (3eV, so that 

n\ = n,2 = fi , (2.66) 

where use has been made of Eqs. (|2.43fl and 1)2.48(1 . 

In order to calculate the thcrmoballistic chemical 
potential n{x) from Eq. I|2.46l) [or, equivalently, the 
thcrmoballistic equilibrium electron density nix) from 
Eq. (|2.48|) ]. the thermoballistic current J(x) from 
Eq. (|2.52|) . and the thermoballistic density n(x) from 
Eq. (|2.53(l . we first have to solve the integral equations 
(12.281) and l|2.31[) numerically for the functions Xi( x ) 
and xz( x )-i respectively. For convenience, we use in 
the integral equations the probabilities p n {x/l) in their 
one-dimensional form, p n ( x/l) = e- x ' 1 [se e Eq. (10) of 
Ref.|23, so that from Eq. igJsj and $FR% 

U n { X ',x";l) = X " ~ X ' e -|« ,, -x , |/t e -e[a a -nun(» / ,»")] 

(2.67) 

(this simplification has only minor effect, see Ref. M)- 
Owing to the scaling properties of the function 
C(x',x";x) given by Eq. (|2.65|) and of the function 
u n (x', x";l) given by Eq. (|2.67() . the results of the calcu- 
lations can be expressed essentially in terms of three di- 
mensionless quantities, for example, x/S, eS, and l/S. In 
Figs. |21 El and^ we show the dependence of n{x), n(x), 



1.5 




FIG. 2: The functions /3[fi(x) — fii] and n(x)/ni, plotted ver- 
sus x/S for eS — 1 and the indicated values of l/S. The ther- 
moballistic chemical potential ^(x) is given by Eq. l|2~4!)jl, and 
the thermoballistic equilibrium density n(x) by Eq. 12.481 . 
The vertical order of the listed l/S values corresponds to the 
order of the curves within the lower and upper set of curves, 
respectively, at x/S = 1. 

J(x), and n(x) on x/S (assuming x\ — 0) for eS — 1 and 
various values of l/S . 

From Fig. |21 (lower panel), it is seen that in the diffu- 
sive limit l/S — > the chemical potential fi(x) decreases 
linearly with x/S and is continuous at the interfaces; 
this behavior persists in the case of arbitrary eS, where 
/3[/J,(x) — fii] = —ex for < x/S < 1. In combination 
with the identical decrease of the potential profile E c (x), 
this implies a position-independent equilibrium density 
n(x) (see upper panel of Fig. [2]). When l/S rises towards 
the ballistic limit l/S — * oo, discontinuities of fi(x) de- 
velop at the interfaces, which increase in magnitude, and 
the slope of /i(cc) becomes smaller. This results in a rise 
of the equilibrium density n(x) across the sample. As a 
function of eS, the discontinuities of /i(x) become smaller 
in magnitude if eS — > oo, and larger if eS — > 0, such 
that in the latter case n(x) becomes independent of x for 
0<x/S<l. 

The ratio J(x)/J shown in Fig. is close to unity 
across the whole sample. As the ballistic limit l/S — > oo 
is approached, J(x) / J becomes more and more symmet- 
ric about x/S — 1/2, and J{x)/J — > 1 in the full range 
< x/S < 1. A somewhat peculiar behavior of J{x) is 
observed in the diffusive regime l/S <C 1. Here, J(x)/J 
is very close to unity inside the sample, except for the 
immediate vicinity of the ends at x/S = and x/S = 1, 
where some structure develops, and J(x)/J converges 
towards a value smaller than unity when x/S — » or 
x/S — > 1. It is important to keep in mind the lat- 
ter feature when, within the unified treatment of spin- 
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1.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x/S 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x/S 



FIG. 3: The ratio J(x)/J, plotted versus x/S for eS = 1 and 
the indicated values of l/S. The thermoballistic current J(x) 
is given by Eq. 112.520 . 

polarized transport, the injected spin polarization at fer- 
romagnet/semiconductor interfaces is defined (see Sec. 
IV. C). When eS is varied, the qualitative behavior of 
J(x)/ J persists for all values of l/S considered. 

In Fig. 0] we show the ratio n(x)/no, where uq is 
the constant value that n(x) [as well as n(x)\ assume 
in the diffusive limit l/S — > 0. For increasing l/S, i.e., 
as the ballistic contribution to the transport mechanism 
increases, the ratio n(x)/no decreases as a whole. This 
effect can be interpreted as reflecting the fact that bal- 
listically, the density decreases rapidly as the velocity 
rises along the sample [see Eq. (|2.1(Jll ]: this is impeded 
at the equilibration points, which lie very dense when 
l/S <C 1 [slow decrease of n(x)] and are widely spread 
when l/S ^ 1 [rapid decrease of n(x)]. For l/S > 1, 
the behavior of n(x) is largely determined by the func- 
tion C(ex) [in the ballistic limit l/S — > oo, we have 
n(x) = (n/2) (1 + e~ eS ) C(ex)]. Again, the qualitative 
behavior of n(x) does not change when eS is varied. 

The zero-bias limit e — * can be treated analytically. 
The solution of Eq. (|2.28|) is found to be 



Xi(x) = 1 + 



x — Xi 
21 



(2.68) 



Since Xi{ x ) = Xi( a; o — x) and Si = ai = 1/2 to zeroth 
order in eS = (3eV, we obtain 



S x-x /2 

x = 1 + 2i' x - [x) = ^zT- 



(2.69) 



From Eq. (|2.52|) . the thermoballistic current J(x) is 
found to be constant, J(x) = J. For the current J, we 
have from Eqs. I2.42|l and l|2.48|l . to first order in (3eV, 

21 

J= 2TTs Veflf3eV > (2J0) 



FIG. 4: The ratio n(x)/no, plotted versus x/S for eS — 1 and 
the indicated values of l/S. The thermoballistic density n(x) 
is given by Eq. fUZfy . and no is the constant value that n(x) 
assumes in the diffusive limit l/S —* 0. 

where n is the common value of the thermoballistic 
equilibrium density at either end of the sample [see 
Eq. H2.66fl ]. Then, the conductance G = eJ/V becomes 



G = 



21 



21 



G, 



(2.71) 



where Q is the Sharvin interface conductance [see Eq. 
(|2.51|) ]. Relation l|2.71|l generalizes the Ohm conduc- 
tance G = (2l/S)G — a/S (valid in the diffusive regime, 
l/S 1), where a = 2f3e 2 v e nl is the conductivity^ The 
thermoballistic equilibrium density n(x) and the ther- 
moballistic density n(x) are obtained from expressions 
(12.48(1 and l|2.53ll . respectively, as 



n(x) — n(x) = n 



(2.72) 



i.e., they are, in the present case of zero bias, both inde- 
pendent of position and equal to the equilibrium density 
at the ends of the sample. 



III. SPIN-POLARIZED TRANSPORT WITHIN 
THE UNIFIED DESCRIPTION 

Having established, in the preceding section, a uni- 
fied description of spinless electron transport in semi- 
conductors in terms of a unique thermoballistic chemical 
potential, we will now extend this scheme by including 
the spin degree of freedom. We allow spin relaxation to 
take place during the motion of the electrons across the 
ballistic intervals. Spin relaxation is generally governed 
by the equation of balance connecting the spin-polarized 
current with the spin-polarized density. In the unified 
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description, it is the thermoballistic current and density 
which enter into this equation. The solution of the bal- 
ance equation is found in terms of a spin transport func- 
tion that is related to the spin-resolved thermoballistic 
chemical potentials. 



A. Balance equation and transport mechanism 

In a stationary situation, the total electron current J — 
Jf (x) + J i (x) composed of its spin-resolved parts Jj (x) 
and J[(x) is conserved, whereas the spin-polarized cur- 
rent J_(x) — J|(x) — Ji(x), or rather its off-equilibrium 
part J-(x) = J-(x) — J-(x), where J-(x) is the relaxed 
part of the spin-polarized current, is connected with the 
off-equilibrium spin-polarized density n_ (x) through the 
balance equation 



d J- (x) fi- (x) 
dx t, 



(3.1) 



Here, n_(x) is defined in analogy to J_(x), and r s is the 
spin relaxation time. For a complete description of spin- 
polarized transport, the balance equation (|3.1I) is to be 
supplemented with a relation between the current and the 
density, which reflects the specific transport mechanism. 

In the ballistic limit, the electron currents Jn(x) are 
proportional to the densities n^i(x) of the electrons par- 
ticipating in the transport, 



Jn( x ) = v ( x ) n n( x ) > 



(3.2) 



where v(x) is the average velocity of the electrons at po- 
sition x [we disregard spin splitting of the conduction 
band edge potential E c (x), so that v(x) is independent 
of spin]. This relation holds also for the off-equilibrium 
spin-polarized current and density, 



J_(x) — v(x)n_(x) . 
Use of this equation in Eq. I|3.1|l yields 



dJ-(x) 
dx 



C{x) 



Is 



o 



(3.3) 



(3.4) 



where l s = 2v e r s is the (ballistic) spin relaxation length, 
which comprises the overall effect of the various under- 
lying microscopic spin relaxation mechanisms ^^S^l^ 
and where Eq. I|2.1U[1 has been used (omitting the posi- 
tions x' , x" of the end points) to express v(x) in terms of 
C(x), i.e., of the potential E c (x) in which the electrons 
move. 

In the diffusive regime, the off-equilibrium spin- 
polarized current and density are connected by the re- 
lation 



J-{x) = — 
e 



dE c {x) 
dx 



1 dfi-(x) 



dx 



(3.5) 



[see Eq. (|2~55|) ]. We then find from Eq. (OJ 
d 2 h-(x) dE c {x) dh-{x) 



dx 2 



dx dx 
d 2 EJx 



where 



dx 2 



n_ (x) 



L 2 



h- (x) = 



(3.6) 



(3.7) 



is the spin diffusion length. 

In the unified description, the total current and den- 
sity inside the semiconducting sample are taken to be 
the thermoballistic current J(x) and density n(x), be- 
tween which no direct relation generally exists. Instead, 
Eqs. (|2.16|) and l|2.17[) express these two quantities sepa- 
rately in terms of the chemical potential fJi(x). The con- 
nection between the off-equilibrium thermoballistic spin- 
polarized current J_ (x) and density n_ (x) can be estab- 
lished along similar lines, as will be described in the fol- 
lowing. 



B. Thermoballistic spin-polarized current and 
density 

In order to include spin relaxation in the unified de- 
scription, we begin by introducing the thermoballistic 
equilibrium densities fi^(x') for spin-up and spin-down 
electrons at an equilibration point x'. It is convenient 
to express fi^(x') in terms of the spin- independent ther- 
moballistic equilibrium density n(x') and a "spin frac- 
tion" a-fj(x') via 



n u( x> ) = n ( x ') a n( x> ) > 



(3.8) 



with aj(x') + ai(x') = 1. In analogy to Eq. i|2.48|l . we 
define spin-resolved thermoballistic chemical potentials 
H U (x') via 



N c e 



-P[E c (x')-n n {x')\ 



(3.9) 



nn( x ') 
which implies 

e /3 m (*') = e W^ atJ (j/) 

The spin fraction a^(x') also enters into the definition 
of the spin-resolved ballistic current Jl,(x',x") injected 
at the left end at x' of the interval \x' , x"], 



(3.10) 



J^(x', x") 



= J l (x',x")a u (x'). (3.11) 

We emphasize that expression (|3.11|) for jl, (x 1 , x") holds 
only at the left end, since this current is not con- 
served owing to spin relaxation, and becomes position- 
dependent inside the interval. There, we write it in the 
form 



J l n {x',x";x) = J\x' ,x")a l u {x' ,x";x) 



(3.12) 
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The function a^(x',x";x) is the spin fraction at po- 
sition x of spin-up (spin-down) electrons injected into 
the ballistic interval [a/, a;"] at its left end at x', with 
a l ^(x' ,x";x) + a l ^(x' , x" ; x) = 1. Here, the dependence 
on the position x" (at the end of the ballistic interval 
[a/, x"] opposite to that at position x' where the electrons 
are injected) is due to the effect of the potential barrier 
embodied in the function E" l (x' , x") in expression l|2.2|l . 

When x coincides with the injection point x' , the cur- 
rent l|3.12[l becomes identical to the current 1)3. 11[) . so 
that 



a\Ax' ,x";x') = a n (x') 



(3.13) 



We now introduce the "spin fraction excess" ot-(x') = 
a-\{x') — ai(x') and the off-equilibrium spin fraction ex- 
cess a-(x') = a-(x') — a~, where &- = a-f — aj., and 
<Sjj are the relaxed parts of the spin fractions (a_ = for 
nonmagnetic semiconductors). With the off-equilibrium 
spin fraction excess a l _ (x' , x"; x) defined in an analogous 
way, we write the off-equilibrium ballistic spin-polarized 
current J l _(x' ,x";x) as 



J l _{x',x";x) = J l (x',x")a l _(x',x";x) . (3.14) 

The spin relaxation of the electrons injected at the left 
equilibration point x' into the ballistic interval [a/, a:"] is 
governed by Eq. I|3.4|) . which, owing to Eq. (|3.12|) . be- 
comes a differential equation for a l _(x' , x"; x), 

da 1 (x',x";x) „, , ., , a 1 (x',x";x) 

- + C(x',x";x)^-^ ^-=0. (3.15) 



dx 

The solution of Eq. (|3~T5l) is 



a l _{x',x"-x)=a-(x')e- lt{x '' x "> x '> x V l ° , (3.16) 



where 



€(x',x"; zi,z 2 ) 



dzC(x',x";z) , (3.17) 



with z < = min(z!,z 2 ) and z > = max^,^)- We then 
have for the off-equilibrium ballistic spin-polarized cur- 
rent at position x of electrons injected at x' 

Jl(x',x";x) = J l (x',x")a.(x')e-^ x '' x "' x '' x ^ ; 

(3.18) 

analogously, we find 

JL(x',x";x) = J r (x',x")a.(x")e^ x '' x "' x - X "^ h 

(3.19) 

for the off-equilibrium ballistic spin-polarized current in- 
jected at x". 

Separating out the relaxed part, we now write the (net) 
ballistic spin-polarized current, in analogy to Eq. 12. 111 , 
in the form 

J-{x',x";x) = J-(x',x";x) + J(x',x")a- , (3.20) 



with the off-equilibrium ballistic spin-polarized current 

J-(x',x";x)=v e N c e-^ x '' x "^ 
x \A (x') e -e(»W.«<)/*. _ A(x") e -e(*'.*";*.*")/».l . 

(3.21) 

here, we have introduced the "spin transport function" 



A(x>) 



-f3[E°-n(x')] 



&_(«') 



(3.22) 



at the equilibration point x' {x\ < x' < x^). Using the 
relation 



1 (\ + a-{rf) 



(3 \l-a^{x') 



(3.23) 



between the splitting (J,-(x') — ^{x') — ^i(x') of the 
spin-up and spin-down chemical potentials and the spin 
fraction excess a-(x'), which follows from Eq. I|3.10|l . we 
have 



A(x) 



-f)[E°-»(x)] 



tanh 



Pfi- (x) 



tanh 



(3.24) 



For the ballistic spin-polarized density, we have, in anal- 
ogy to Eq. JHI231, 

n._(a/, x"; x) — h-(x', x"\ x) + n(x', x"; x) <S_ , (3.25) 

where 

h-(x',x";x) = ^C(x',x";x)e-^( x '' x " ) 

X A{x') e -C(x',s";x',x)/i s + A(x") e -£W,*"w")/l* 

(3.26) 

is the off-equilibrium ballistic spin-polarized density [see 
Eqs. (|2~TT)) and l|3~2TJ) ]. 

For the thermoballistic spin-polarized current J-(x) 
passing through the point x, we find 



J_ (x) ~ J_ (x) + J(x) a_ , 



(3.27) 



where the off-equilibrium thermoballistic spin-polarized 
current J_(x) is obtained from the off-equilibrium bal- 
listic spin-polarized current l|3.21|l by summing up the 
weighted contributions of the ballistic intervals, 



J-(x) = v e N c 



dx' 

T 



dx' 1 

T 



(x',x";l) 



A(x') e~^ x ' > x "\ x> _ A(x") e~ e ( x ',x"-,x,x")/l s 

(3.28) 

(a~! < x < x 2 ). Similarly, the thermoballistic spin- 
polarized density n_(x) at the point x is 



n_ (x) — fi- (x) + n{x) <S_ , 



(3.29) 



12 



where the off-equilibrium thermoballistic spin-polarized 
density n_ (x) is obtained from l|3.26[) as 



, , N c [ x dx 1 f X2 dx" , , .. , x 
h_(x) = ^j —J —Wc(x',x";l;x) 

A{x') e -e(.x',x"-y,x)/t, + A ( x »} e -€(x',x";x,x")/l s 

(3.30) 

(xi < x < x 2 )- 

In the diffusive regime, l/l s <C 1, l/S -C 1, the integrals 
over x' and x" in Eqs. (j3.28(l and (|3.30|l can be evaluated 
explicitly, yielding 



J_(ar) =-2v e N c le'^ E ^- E ^ 



dA(x) 
dx 



(3.31) 



and 



= 7V C e -/S[«.(»)-B?] A(ar) . (3.32) 



Eliminating the function A(x) from these two equations, 
we obtain the standard drift-diffusion relation l|3.5|l . 

In the current J-(x) and density h—(x), the spin re- 
laxation in each ballistic interval is described in terms of 
the values of the spin transport function A(x) at the end 
points x = x' and x = x" . Since J-{x) and fi-(x) are 
linearly connected with A(x), they are linearly connected 
with each other. Thus, it appears that the function A(x) 
[and not the chemical-potential splitting is the 

key quantity for treating spin transport within the uni- 
fied description. It remains to find an equation for the 
determination of this function. 



C. Integral equation for the spin transport function 

The required equation is provided by the basic bal- 
ance equation, Eq. I|3.1|l . which we now read in terms of 
the off-equilibrium thermoballistic spin-polarized current 
(|3.28l) and density l|3.30|l of the unified transport descrip- 
tion. Since the derivative with respect to x of the terms 
in the brackets of expression l|3.28|l for J_ (x) is compen- 
sated by the term n_(x)/r s (this reflects the fact that 
spin relaxation in the ballistic intervals has already been 
taken into account), only the derivative on the limits of 
integration in expression i|3.28f) remains, and we have 



X2 dr' r 

— W(x,x';l) A(x) - A{x') e- €( - x ' x ^ u 



dr' v 

— W{x', x- 1) \ A{x') e-^ x ■*>/'- - A(x) 



where 



£{x',x") = £(x',x";x',x") 



= 0, 
(3.33) 

(3.34) 



With the action of the symbolic operator W(x' ,x";l) 
explained by comparison of Eqs. (|2.12|) and 12.15JI . 
Eq. H3.33J) reads explicitly 

W 2 (xi,x- 1 l,l s )A 1 + W 2 (x,x 2 ;l,l s )A 2 
- W(x!,x 2 ;x;l) A(x) 
f X2 dx' 

+ / ^W 3 (x' ) x;l,l s )A(x')=Q, 

(3.35) 

where 

W n (x', x"; I, l s ) = w n (x', x"; I) e ^ x '- x "^ 1 ^ , (3.36) 



W(x\,x 2 ;x; I) 



w 2 (xi,x;l) 
X2 dx[ 
, I 



w 2 (x,x 2 ;l) 
w 3 (x',x; I) , (3.37) 



and Ai.2 = A(xi y2 )- Equation (|3.35l) is a linear, 
Fredholm-type integral equation for the spin transport 
function A{x). Its solution for x\ < x < x 2 determines 
the spin-polarized electron transport inside the semicon- 
ducting sample, and is obtained in terms of the values 
A\ and A 2 at the interfaces at the ends of the sam- 
ple. The latter are determined by the off-equilibrium 
spin fraction excesses a-(x\^) and the chemical poten- 
tials fj,\ t 2 = n(xi y2 ) in the contacts [see Eq. (|3.22|l ]. The 
function A{x) is not, in general, continuous at the inter- 
faces, A(xi) 7^ Ax, A(x^) 7^ A 2 ("Sharvin effect"), as 
will be demonstrated in Sec. III.D by way of a particu- 
lar example. The discontinuities of A(x) arise from the 
joint effect of the discontinuities of the spin-independent 
thermoballistic chemical potential /i(x) and those of the 
spin fraction excess a-(x) [or, equivalently, of the spin- 
resolved thermoballistic chemical potentials /i-\i(x)]. 

Substituting A(x) in Eqs. lj3~2~5)) and ijOty . we obtain 
the off-equilibrium thermoballistic spin-polarized current 
J~(x) and density n_(x), respectively; the thermoballis- 
tic spin-polarized current J_ (x) and density n_ (x) then 
follow from Eqs. I|3.27ll and 13.291) . respectively. Dividing 
by the corresponding total thermoballistic current and 
density, Eqs. (|2.16|) and (|2.17|) . respectively, we get the 
current spin polarization 



PA*) 



J-(x) _ J-(x) 
J(x) ~ J{x) 



and the density spin polarization 

Tl-(x) fl-(x) 
Pn{X) = = + (X- 

n(x) n(x) 



(3.38) 



(3.39) 



inside the sample. These polarizations are written in 
terms of the thermoballistic current and density; how- 
ever, we take their magnitudes to be also those of the 
physical polarizations, for the following reason. The un- 
derlying assumption of our approach is that the equilibra- 
tion process, i.e., the coupling between the thermoballis- 
tic and background currents, is independent of spin (this 
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is clearly true for the D'yakonov-Perel' spin relaxation 
mechanism^i but remains to be examined for the other 
mechanisms). Therefore, the relative spin content is the 
same in these two currents, and thus equal to that of their 
sum, viz., the physical current. Hence, wc may take the 
polarizations Pj(x) and P n (x) of Eqs. H3.38J) and H3.39J) . 
respectively, for the physical polarizations. 

The integral equation l|3.35|l constitutes the central re- 
sult of the present work. It allows the calculation of the 
spin polarization in semiconductors for any value of the 
momentum and spin relaxation lengths as well as for ar- 
bitrary band edge potential profile. The fact that we are 
led, in the unified description of spin-polarized transport, 
to an integral equation is connected with the introduc- 
tion of the momentum relaxation length / as an inde- 
pendent parameter of arbitrary magnitude, which gives 
rise to nonlocal ballistic effects. The basic parameters 
controlling the transport in the unified description are 
the equilibrium densities the momentum relaxation 
length I, and the spin relaxation length l s , whereas in the 
standard drift-diffusion model one uses the conductivity 
a and the spin diffusion length L s . 



D. Differential equation for the spin transport 
function 



In order to interpret our unified description of spin- 
polarizcd transport and relate it to previous, less gen- 
eral descriptions, we consider in the following the case of 
field-driven transport in a homogeneous semiconductor 
without space charge. As in Sec. II. C, we take the prob- 
abilities p n (x/l) in their one-dimensional form. Then, 
Eq. H3.35fl can be converted into an integrodifferential 
equation for the spin transport function A(x). In an ap- 
proximation which is adequate for the present purposes, 
the latter equation reduces to a second-order differential 
equation. 



1. General form and diffusive regime 

For a potential of the form l|2.64[l , the integral equation 
(THB) redu ces, with the help of Eqs. (ESSjl, 1233, GUI, 
and (jOHl, to 

fi(x -Xi)Ai+ f 2 (x 2 - x) A 2 

r x dx' 

- f(x-xi)A(x)+ / — fl ( x - x >) A(x') 



where 



f X2 dx' 

J ^-f 2 (x'-x)A(x')=Q, 

fax) = e -N-i/J+e(«)/J.]* . 
f 2 {x) = e -Wl+^)/ls]x } 



(3.40) 



(3.41) 



(3.42) 



/(*) = 



1 + d 



-(e+l/l)x 



and 



C(C) = C7o d( ' CiC) 



(3.43) 



(3.44) 



with < c(C) < 1, c(C) -> 1 for C -> 0, and c(C) ~ 
2(?<r 1/2 for C -> oo. 

By supplementing the inhomogeneous integral equa- 
tion (|3.40[1 with the equations obtained by forming its 
first and second derivative with respect to x, and elimi- 
nating from this set of equations the quantities A\ and 
A 2 , we can convert Eq. I|3.40|) into a homogeneous in- 
tegrodifferential equation for the spin transport function 
A{x). [This procedure could also be applied to Eq. H3.35J1 . 
but does not seem to be helpful in the general case] . Now, 
the latter equation can be simplified by replacing the 
function c(£) with a position-independent average value 
c, so that the coefficient functions fi(x) and f 2 (x) in 
Eq. (|3.4l)|l reduce to pure exponentials. With this ap- 
proximation, the integrodifferential equation for A(x) be- 
comes a second-order differential equation of the form 



, , . d 2 A{x) , , . dA{x) 



dx 2 



dx 



+ b 2 (x)A(x) =0, (3.45) 



where 



6o(af) = 2 + el[l + b(x)} , 



b 1 {x)=e(2 + el)[l-b(x)] 



(3.46) 



(3.47) 



2 (x) = -Ll 



i 2 -u + di{i - 1) 



el(l-l)(l + l + ell)[l + b(x)]} , (3.48) 



with 



and 



b(x) = e -(«+i/0(«-»0 



1 1 



I I, 



1 7 _ Is 
~ j 's — - ■ 
C 



(3.49) 



(3.50) 



Since, owing to the presence of the factor e~ x ^ L in the 
functions fi(x) and f 2 (x), only the values of c(ex) within 
the range < x ^ I contribute appreciably, we choose c 
as the average of c(ex) over an x-interval of length equal 
to the momentum relaxation length I, 



I 



dx c(ex) 



d 



d( \n(d/() C(C) • (3.51) 



In the right-hand integral of this equation, the range of 
small ex, where C(ex) sa 1, is emphasized because of the 
weight factor ln(d/£). For large el (in the ballistic regime 
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and/or for strong fields), the variation of c(ex) with x be- 
comes essential, and a more detailed study of the validity 
of the approximation leading to Eq. 1(3.45|1 will be neces- 
sary. For the present purpose of solely demonstrating the 
principal effects of the transport mechanism, we consider 
this approximation, in conjunction with the choice (|3.51|1 
for c, to be sufficiently accurate. 

In the diffusive regime characterized by the conditions 
l/l s <C 1, l/S <C 1, and el <§; 1, we have c = 1 and I = I, 
and Eq. (|3.45l) reduces to 



(3.52) 



1 



d 2 A(x) 



dA(x) 1 
dx 2 • £ —-Ll A{x)=0 



In view of Eq. I|3.32|) . Eq. (|3.52l) can be rewritten in 
terms of n_(x) and then agrees with Eq. I|3.6[) . and with 
Eq. (2.8) of Yu and Flatted if the intrinsic spin diffu- 
sion length L of that reference is identified with the spin 
diffusion length L s = yllg. It thus turns out that Eq. 
(13.45(1 generalizes the usual spin drift-diffusion equation 
to the case of arbitrary values of the ratio l/l s . 



2. Zero-bias limit 



In the zero-bias limit e 
(|3.40j) reduces to 



0, the integral equation 



e -{x-X-L)/lA _|_ e -(x 2 -x)/l £ 

- 2A(x) + / ^L e -\ x ~ x \' l A{x') 



, (3.53) 



from which one derives the differential equation 



d 2 A(x) 
dx 2 



here, 



L 



y/l + l/h 



(3.54) 



(3.55) 



is the generalization of the spin diffusion length l|3.7|l . 
which includes ballistic effects via the renormalization 
factor 1 / y/l + l/l s . The length L becomes equal to the 
spin diffusion length proper, L = L s , in the diffusive 
regime where 1 = 1, and to the spin relaxation length, 
L = l s , in the ballistic limit ljl s — > oo where I = l s . 

Equation (|3.54|) has the general solution, for 27 < x < 
x 2l 



A(x) = C ie - {x - Xl)/L + C 2 e~ {x2 ' x)/L 



(3.56) 



With this expression substituted for A(x) in Eq. I|3.53|l . 
the set of two equations resulting from writing down this 
equation for x = x% and x = X2, respectively, can be 
solved for the coefficients C\ 2, 



1 r 
D 



(l+ 7 )e s/L A 1 -(l- 7 )A 2 



(3.57) 



D 



(l- 7 )^i-(l+7)e S/L A 2 



(3.58) 



C 2 = 
where 

D= (l+ 7 ) 2 e s / L - (l- 7 ) 2 e- s / L , (3.59) 



with 



7 



L I 



L L V l+l 



< 1 . 



(3.60) 



It follows from Eqs. (|3.56|l - (|3.60l) that the function A(x) 
is discontinuous at x = x\^, 

AAi = A(x+) - A x = ~ { 9 A 1 - hA 2 ) , (3.61) 



1 



-(hA 1 -gA 2 ) 



(3.62) 



AA 2 = A 2 - A(x 2 ) 
where 

[(1 + 7) e s/L + (1 - 7) e- s ' L ] < 1 (3.63) 



27 



g= D 



and 



D ~ I + 7 



(3.64) 



In the diffusive regime, one has L = L s = l s -Jl/l s and 
7 = 0, and therefore 

A{x) = Ax e - [x ~ x ^l Ls + A 2 e -i^-x)/L 3 ^3 65 ) 
In the ballistic limit, one has L = l s and 7 = 1, so that 



A{x) 



1 r 

2 



Ai e -( X - X1 )/l s + ^ e -(x 2 -x)/l s 



(3.66) 



The discontinuity of A(x), e.g., at x = xi, is AAi = 
A 2 exp(—S/L s ) in the diffusive regime, and AA% = 
h[— At + A 2 exp(— S/l s )] in the ballistic limit. 



IV. SPIN-POLARIZED TRANSPORT IN 
FERROMAGNET/SEMICONDUCTOR 
HETEROSTRUCTURES 

We now turn to the unified description of spin- 
polarized electron transport in heterostructures formed of 
a semiconductor and two ferromagnetic contacts (cf. Fig. 
0. We treat the ferromagnets as fully degenerate Fermi 
systems. The semiconductor is taken to be nonmagnetic 
[i.e., <5_ = 0, and hence a^(x') = a-(x')] and homoge- 
neous without space charge. We disregard spin-flip scat- 
tering at the interfaces, but spin-selective interface resis- 
tances are included in our description by introducing dis- 
continuities into the spin-resolved chemical potentials, in 
the same way as in previous descriptionsi^iSiiiiS within 
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the drift-diffusion model. Of course, for realistic applica- 
tions, it is necessary to treat the effect of interface barri- 
ers explicitly^ using potential profiles E c (x) which, gen- 
erally, must be calculated self-consistently from a nonlin- 
ear Poisson equation. In some cases, however, it may be 
sufficient to perform non-self-consistent calculations us- 
ing appropriately modeled band edge profiles In 
any case, this would require the spin transport function 
A{x) to be determined by numerically solving the inte- 
gral equation (|3.35|) . This task will be deferred to future 
work. 

In order to obtain the position dependence of the spin 
polarization across the heterostructure, the current spin 
polarization and the chemical potential in the semicon- 
ductor are to be connected with the corresponding quan- 
tities in the left and right ferromagnets. The current spin 
polarization Pj(x) [see Eq. (|3.38(l ]. as expressed by the 
ratio of the thermoballistic currents J_(x) and J(x), is 
equal to the physical current spin polarization [see the 
discussion following Eq. (|3.39() ]. It is, therefore, contin- 
uous across the whole heterostructure, in particular, at 
the interfaces, so that the current spin polarizations in 
the semiconductor and the ferromagnets can be equated 
directly there. On the other hand, in the presence of 
ballistic contributions, the thermoballistic chemical po- 
tential //(x) and the spin transport function A(x) are 
not continuous at the interfaces. The discontinuities at 
the interfaces are taken into account when the functions 
and A(x) inside the semiconductor are calculated 
in terms of their values and A\ %2l respectively. The 
latter values are to be equated with the values of the 
corresponding quantities in the ferromagnet. 

We begin with a brief summary of the standard de- 
scription (see, e.g., Ref. IT^) of the spin polarization in 
the ferromagnets. 



A. Current spin polarization in the ferromagnets 

In the (semi-infinite) left ferromagnet located in the 
range x < X\, the spin-up and spin-down chemical po- 
tentials ^ii(x) are given by 



Ml ( x ) 



0"! 



r (1) 



(4.1) 



we now find for the current spin polarization 
Gi 



Pj(x) =P 1 - 



where Pi — (cr j 1 ^ — o~Y' > )/o~i is the relaxed (current or 
density) spin polarization in the left ferromagnet, and 



(4.3) 



G x = 



4^f> 



0\ L 



(4.4) 



is a transport parameter of the ferromagnet, which has 
the dimension of interface conductance. Analogously, we 
obtain 

Pj{x) = p -2 + ^Fj»- (4) e-^^ (4.5) 

for the current spin polarization in the right ferromagnet 
located in the range x > x 2 . 

In the absence of spin-selective interface resistances, 
the chemical-potential splitting /i_(x) is continuous at 
the interface, /i~(x^) = /j,-(xi) and fi^(x 2 ) — /j.-(x 2 ), 
where fi-(xi t2 ) are its values at the interface itself. The 
latter are to be set equal to the corresponding values in 
the semiconductor, which yields 

[M— (^l,2)]ferromagnct — [M— (*£l,2)]scmicoiiductor 

-3"»(1±S54), (4.6) 

p \^-ai,2J 
where the right-hand part of this equation follows from 
Eq. I|3.23[) for %' = xi,2, and «i,2 = OL— (^1,2)- F° r the 
current spin polarizations at the interfaces, Pj(xi^), we 
have from Eqs. I|4.3|) and l|4.5(l 



Pj{xi) = Pi - „ „ In 



Pj{x 2 ) =P 2 + 



2/3e 2 J 

G 2 
2/3e 2 J 



In 



1 — a?i 



1 + a 2 

1 - a 2 



(4.7) 



(4.8) 



which are to be set equal to the corresponding polariza- 
tions of the semiconductor. 

(1 2) 

Spin-selective interface resistances pi, 1 are introduced 
via discontinuities of the spin-resolved chemical poten- 
tials on the contact sides of the interfaces. At x = xi, 
for example, the discontinuity has the form 



P (1) 

ni 



(4.9) 



The corresponding interface resistance is located between 
x = x^ and x\ (in the ferromagnetic contact), and thus 
where Lf is the spin diffusion length. The quantities i s adjacent to the Sharvin interface resistance between 

x\ and x\ (in the semiconductor). The quantity 



crj^ are the conductivities for spin up and spin down 



which are independent of position, and a± = ai + a. 
We then have C\ = ai 1 ' /if (2^) = — a^fj,i(x^), and 
therefore C\ = (a^a^/ai) /^(x^), where p~(x) — 
M( x ) - Ml (a)- With 



(i) 



fj.-(x 1 ) to be substituted in Eq. (|4.3() is obtained, using 
Eqs. lpTT)l . g2J>, and ljO|) . as 



fj,_(x 1 



1 



GrpW/4 



Jn( x ) 



e 2 dx 



(4.2) 



2 J 



-(*!) + — 



P 



(4.10) 
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where = ± p^ 1 ■ The connection of p-(xi) with 
the spin fraction excess a\ is, as before, given by Eq. 
(|4.6|) . The same procedure applies mutatis mutandis to 
the interface at x = x 2 . 



B. Spin polarization across a heterostructure in 
the zero-bias limit 



In the zero-bias limit J — > 0, we now demonstrate the 
procedure for calculating the current and density spin 
polarizations across a ferromagnet/semiconductor het- 
erostructure. 

Evaluating expressions 1)3. 28[) and (|3.3U|) . respectively, 
with A(x) given by Eq. I|3.56[) . we find for the thermobal- 
listic spin-polarized current in the semiconductor 



J_(x) 



-dA(x) 
dx 



(4.11) 



and for the thermoballistic spin-polarized density 

n_(x) = N c A(x) (4.12) 

[x\ < x < x%). For zero bias, one has J(x) = J = 
const, and n{x) = n = const, [see Eqs. ()2.70(l and l|2.72(l ]. 
so that, by combining expressions (|4.11|l and (|4.12|) . we 
obtain the relation 



PA*) 



2v e nl dP n (x) 



J 



dx 



(4.13) 



between the current and density spin polarizations. Fur- 
thermore, J_(x) and n_(x) both satisfy Eq. (|3 . 54|) . 
and so do the polarizations Pj(x) and P n {x) given by 
Eqs. (|3.38J) and (|3.39(l . respectively. Differentiation of 
Eq. (|4"llH) then yields, together with Eq. IpTo!)) for P n {x), 



Pn(x) = - 



Is J dPj( X ) 



2v P n dx 



(4.14) 



From Eq. I|3.38|l with J(x) = J, we find, using Eqs. (|3.5GI) 
and (|4.11(l , the explicit form of the current spin polariza- 
tion as 



2v e NJ 
LJ 



q e -(x-x t )/L _ q e -{x 2 -x)/L 

(4.15) 



The density spin polarization is obtained from Eqs. 
and H4.12fl [or, equivalently, from Eqs. (|4.14() and 14.15|l ] 

as 



P n {x) 



Nc 
n 



Ci e -( x - Xl)/L + C 2 e ~ {xo -- x)/L 



(4.16) 



The coefficients C\^ 2 in Eqs. (|4.15() and l|4.16[) can be ex- 
pressed via Eqs. H3.57[H3.60|) . using Eq. I|3.22|l . in terms 
of the spin fraction excesses ol\ i2 on the contact sides of 
the interfaces. 



In order to determine the quantities otx,2, we consider 
the current spin polarization l|4.15|l on the semiconductor 
sides of the interfaces, 



Pj(xt) = a o T (d a i - ha 2 ) 



Pj{*2) 



l3e 2 J 



l3e 2 J 



{hai - ga 2 ) ; 



(4.17) 



(4.18) 



here, Q is the Sharvin interface conductance given by 
Eq. (|2.51() . and the coefficients g and h are given by 
Eqs. I|3.63|l and (|3.64(l . respectively. As mentioned be- 
fore, Pj(xi) — Pj{xi) and P.jix^) — Pj(x 2 ), and the 
connection with the polarization in the contacts is made 
by equating expressions (|4.7|l and (|4.17|l . and expres- 
sions (|4.8(l and (|4.18|l [note that, if spin-selective inter- 
face resistances are included, expression (|4.7|) for Pj(xi) 
is to be replaced with the general expression obtained 
by using expression l|4.10|l for /i_(x^) in Eq. (|4.3() . and 
analogously for Pj(x 2 )]. In the zero-bias limit, when 
J(xia) — J [or k = 1 in Eq. (|2.24|l ]. we have |ck 1j2 | < 1, 
and this procedure then results in the system of coupled 
linear equations 



(g + Gi) 



Gi) a\ — ha 2 



Pe 2 J 



Pi , 



hai - [g + G 2 ) a 2 = 



(3e 2 J 



Pi , 



where 



Gi, 



(4.19) 



(4.20) 



(4.21) 



The solutions of Eqs. (|4.19|) and l|4.2U|l are found to be 



Oil 



a 2 



(g + G 2 ) 



Pi - hP 2 



where 



(3e 2 J 



pe 2 J 



A= (g + Gi) ig + G 2 )-h 



hPt-lg + Gx)^ 



(4.22) 



(4.23) 



(4.24) 



Expressions l|4.22|l and l|4.23|l determine the spin frac- 
tion excesses a\ and a 2 in terms of the current J, of the 
polarizations P± and P 2 in the left and right ferromag- 
net, respectively, and of material parameters, such as the 

(1 2) 

conductivities ai i2 and the spin diffusion lengths L s : ' 
of the ferromagnets (via Gi i2 ), and the momentum re- 
laxation length I and the spin relaxation length l s of the 
semiconductor as well as its length S (via g and h) and 
the equilibrium density n (via Q). Since the quantities 
cm 2 are proportional to the current J, the current spin 
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x [nm] 

FIG. 5: The zero-bias (e — > 0) current spin polarization Pj(x) 
along a symmetric ferromagnet/semiconductor/ferromagnet 
heterostructure with S = 1 /im for the indicated values of the 
momentum relaxation length /, calculated from Eqs. 14. 3L 
14.511 . and 114.1511 with X\ = and X2 = S. The solid curves 
correspond to zero interface resistance, the dashed curves to 
interface resistances of 10~ 7 fl cm 2 for spin-up electrons and 
5 x 1CT 7 n cm for spin-down electrons, respectively. For the 
remaining parameter values, see text. 

polarization Pj(x) is independent of J, while the density 
spin polarization P n (x) is proportional to J. 

The current spin polarization along the entire het- 
erostructure, Pj(x), is now obtained as follows. In the 
semiconductor, it is given by expression (|4.15(l . with C12 
calculated from as explained there. In the ferro- 
magnets, the expressions for the current spin polariza- 
tion are provided by Eqs. (|4.3[) and l|4.5l) . respectively, 
where the quantities /i_(xj~) and /i_(x^~) are calculated 
from Eq. (|4.1U|) and from its analogue for //_ (xf), respec- 
tively. Analogously, the density spin polarization P n (x) 
in the semiconductor is given by expression Q4.14)) . We 
do not write down the density spin polarizations in the 
ferromagnets, but only mention that they do not, in gen- 
eral, match the polarizations P n (x±) and P„(x^~) on the 
semiconductor sides of the interfaces. 

In order to demonstrate the effect of the trans- 
port mechanism (characterized by the magnitude of 
the ratios l/l s and l/S), we show in Fig. [5] the zero- 
bias current spin polarization Pj(x) for a symmet- 
ric ferromagnet/semiconductor/ferromagnet heterostruc- 
ture with sample length S = 1 /im at T = 300 K as a 
function of x for various values of the momentum relax- 
ation length I. For the parameters of the ferromagnets, 
we adopt from Ref. the values <j\ = 02 = 10 3 fi" 1 
cm -1 for the bulk conductivities and Li' = L s 2 ^ = 60 
nm for the spin diffusion lengths; the bulk polarizations 
are chosen as P\ = P2 = 0.8. For the material parameters 



of the semiconductor, we take the values m* = 0.067m e 
for the effective electron mass, l s — 2.5 /im for the bal- 
listic spin relaxation length (corresponding to n-doped 
GaAs; see Refs. and El), and n — 5.0 x 10 17 cm~ 3 
for the equilibrium electron density. Clearly, in a specific 
semiconducting system, the value of the momentum re- 
laxation length I is fixed. Therefore, when varying I, we 
are considering the above parameter values to be repre- 
sentative for a whole class of semiconductors (regarded 
as nondegenerate; at room temperature, this should be 
an acceptable working hypothesis^ 3 -) that differ in the 
strength of impurity and phonon scattering and hence in 
the magnitude of I. 

The momentum relaxation length I affects the results 
shown in Fig. [5] in a twofold way. (i) It determines the 
conduction in the semiconductor. For small values of I, 
the conductance of the latter is small, and thus the con- 
ductance mismatch with the ferromagnets is large, lead- 
ing to a small injected current spin polarization Pj(0). 
(ii) It determines the generalized spin diffusion length 
L = [ll s /(l + l/ls)] 1 / 2 , which acts as the polarization de- 
cay length, so that for small I the polarization dies out 
rapidly inside the semiconductor. The degree of polar- 
ization may be raised considerably all along the semicon- 
ductor when the value of I is increased up to a length of 
the order of the sample length, in which case the ballistic 
component becomes prevalent. Figure|S]also shows that, 
by introducing appropriately chosen spin-selective inter- 
face resistances, one may offset the suppression of the 
injected polarization due to the conductance mismatch 
for small I; however, the rapid decay of the polarization 
inside the semiconductor cannot be prevented in this way. 

For the case of Fig- EJ we show in Figs. [5] and respec- 
tively, the zero-bias current spin polarization Pj(x) for 
various values of the equilibrium density n and the spin 
relaxation length l s . It is seen that varying n has about 
the same overall effect on Pj(x) as varying the momen- 
tum relaxation length I, whereas varying l s affects mainly 
the rate of decay of Pj(x). 



C. Injected spin polarization for field-driven 
transport 

We introduce the "injected spin polarization" as the 
spin polarization at one of the interfaces, e.g., at x = x\, 
generated by the bulk polarization P\ of the left fer- 
romagnet regardless of the influence of the right ferro- 
magnet. More precisely, we define the injected current 
spin polarization as the current spin polarization Pj(x±) 
given by Eq. (|3.38[) in the limit S/L — > 00. Similarly, 
the injected density spin polarization is defined as the 
polarization P n (xf) of Eq. (|3.39|l in the same limit. The 
injected spin polarization at x — x^ provides the initial 
value of the left-generated polarization in the semicon- 
ductor, which propagates into the region x > x\ while 
being degraded by the effect of spin relaxation. 

We now consider the injected spin polarization for elec- 
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FIG. 6: The zero-bias current spin polarization Pj(x) for the 
case of Fig. |^| for zero interface resistance and various values 
of the momentum relaxation length / and of the equilibrium 
density n. Short-dashed curves: n — 4 x 10 17 cm -3 ; full 
curves: n = 5 x 10 17 cm -3 ; long-dashed curves: n = 8 x 10 17 

™-3 



tron transport driven by an external electric field, i.e., a 
potential profile of the form 



1. General case 

In order to obtain the spin transport function A(x), 
we have to solve Eq. (|3.45[) numerically under the con- 
dition A(x) (x exp(— x/X) for x — > oo. The decay 
length A is determined by solving Eq. I|3.45[l in the range 
x — x\ ^> (e + l/l)^ 1 where the function b(x) in the 
coefficient functions bo(x), b\{x), and 62 (%) can be disre- 
garded, 



A = < ~ 



1 + el 1 + el 2 + el 



I'- 



ll 2 + el 



1/2' 



(4.25) 



It can be shown that for any combination of parameter 
values, A is a real number. 

For calculating the injected current spin polariza- 
tion from Eq. I|3.38[) . we determine the thermoballis- 
tic spin-polarized current at the interface, J-(xi), from 
Eq. Using Eq. I|3.22[) and fixing the normaliza- 

tion of the function A(x) in terms of A\ with the help of 
Eq. ESS) , we find 



J-{x~l) = v e nT j ai , 



(4.26) 



FIG. 7: The zero-bias current spin polarization Pj{x) for the 
case of Fig. |S| for zero interface resistance and various values 
of the momentum relaxation length I and of the spin relax- 
ation length l s . Short-dashed curves: l 3 = 1 /im; full curves: 
Z s = 2.5 /tm; long-dashed curves: h = 10 /im. 



where 



and 



A(xf) - A 
A(xf) - A/2 



A 



dx 
T 



-(x—xi_)/l 



A(x). 



(4.27) 



(4.28) 



To find the total thermoballistic current at the interface, 
J(xf ), we go back to Eq. (|2.24|l . Expressing the current 
J in the form 



1 

X 



J = — v e n , 



(4.29) 



which follows, for e > and S/L — > 00, from the 
current-voltage characteristic l|2.42|l with (3eV ~ eS and 
N c exp(—[3Ep) = n, we obtain 



J(xi) — — v e n 



(4.30) 



This expression is conveniently evaluated by using for k 
and x the closed-form representations 



1 + d 

2 + el 



X 



(1 + eQ 2 
el{2 + el) 



(4.31) 



which have been inferred from the results of systematic 
numerical calculations for fixed el > and very large 
values of S/L. For the injected current spin polarization, 
we now find 



PA4) 



J(xf) 



X 



rjai , 



(4.32) 
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FIG. 8: The injected current spin polarization Pj(xi) for 
S/L — > oo as a function of the electric-field parameter e for 
the indicated values of the momentum relaxation length I. 
The solid curves correspond to zero interface resistance, the 
dashed curves to interface resistances of 10 -7 Q cm 2 for spin- 
up electrons and 5 x 10 fi cm 2 for spin-down electrons, 
respectively. For the remaining parameter values, see text. 

which, by continuity, is equal to Pj(x\). Setting the 
right-hand side of Eq. I|4.32[l equal to expression 1)4.711 [or 
to the more general expression including spin-selective 
interface resistances; see the remark following Eqs. (14.17(1 
and (|4.18|l ] for the injected spin polarization in terms of 
the contact parameters, we arrive at 

Pr-^Xln(\^)=^T J a 1 . (4.33) 
2 \i — atxj re 

This is a nonlinear equation for ot\ which is to be solved 
for given values of the parameters e, Pi, Gx, n, I, and l s . 

Turning to the calculation of the injected density 
spin polarization, we determine the thermoballistic 
spin-polarized density at the interface, n_(a;^~), from 
Eq. IpHUj) . using again Eqs. (|3~2~2T) and lp£lTj|) . and obtain 



n>-( x t) = - T n ax 



where 



r n = 



A(x+) 



A(x+) - A/2 



(4.34) 



(4.35) 



For the total thermoballistic density at the interface, 
n(xf ), we find from Eq. using Eqs. ij^igjl . lt232l . 

and (Og , 



n(xi) = n I 1 



K.J 

2v P h 



= n 1 

2 X 



(4.36) 



The injected density spin polarization now follows as 



Pn{xt) 



n(xj) 2\ — K 



T„ ax 



(4.37) 



where the spin fraction excess ax is again to be deter- 
mined by solving Eq. I|4.33|l . 

Figure |H1 shows the injected current spin polarization 
Pj{xx) for S/L — ► oo as a function of the electric-field 
parameter e for various values of the momentum relax- 
ation length I; the remaining parameter values are the 
same as in Fig. El In calculating Pj(xx) from Eq. I|4.32|l . 
we have used k = 1. This choice has been made because 
expression (|4.32() with k given by Eq. I|4.31|l does not 
represent a meaningful injected polarization in the diffu- 
sive limit (see below) ; instead, one must set n = 1 in this 
limit. For simplicity, we have used this value throughout. 
In conformity with the drift-diffusion results of Ref. figt 
the injected polarization generally rises with increasing 
e; however, as in Fig. \E\ the main effect is due to the 
variation of I. 



2. Diffusive regime 

In order to relate our treatment of the injected spin 
polarization at ferromagnet/semiconductor interfaces to 
previous treatments within the drift-diffusion model, 
in particular to that of Yu and Flatted we consider 
the diffusive regime, l/l s <C 1 and el -C 1, in some 
detail. In that regime, the spin transport function 
A(x) is determined by Eq. (|3.52|) . whose solution is 
A(x) oc exp(— x/Ll), with the field-dependent spin diffu- 
sion length L\ given by 



1 



1 



1/2 



(4.38) 



We then obtain Tj = 21 /L\ and T n = 2. Furthermore, 
from Eqs. I|4.31|l for el <C 1, we find k = 1/2 and x = 
l/2eZ. 

At this point, some analysis is required regarding the 
definition of the injected current spin polarization in the 
diffusive regime. In the definition introduced above, 
first the functions J_(x) and J(x) are evaluated for 
x — * Xi, and subsequently the diffusive limit is ap- 
proached. This procedure results, in particular, in the 
value k = J(x^)/J —1/2. A closer look at the function 
J{x)/J (see Fig. OJl, however, shows that in the diffu- 
sive regime this function is virtually equal to unity inside 
the semiconducting sample and tends to smaller values 
only within a (very short) distance of order I from the 
interfaces. Therefore, it is indicated here to define the 
injected current spin polarization in terms of a position 
x > I inside the sample, where J(x)/J — 1 is the rel- 
evant value for the propagation of the spin polarization 
into the semiconductor. Thus, in the diffusive regime, 
we adopt the effective value re = 1 in the calculation of 
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the injected spin polarization. While in the ballistic limit 
(and now also in the diffusive regime) the choice re = 1 
is unique, in the range of intermediate Z-values a mean- 
ingful definition of the injected spin polarization requires 
an appropriate choice of the position inside the sample 
at which the thermoballistic current and spin-polarized 
current are to be evaluated. 

With the choice re = 1, the injected current spin polar- 
ization in the diffusive regime is obtained from Eq. (|4.32() 
as 



(4.39) 



where the spin fraction excess ol\ is now to be calculated 
from Eq. (|4.33|) with re = 1. Since re < x I0r d < 1; tne 
injected density spin polarization in the diffusive regime 
follows from Eq. (|4.37|l as 



Pn (^1 



Oil 



(4.40) 



where ai again must be calculated from Eq. H4.33[l with 
re = 1. 

Comparing our results for the injected spin polariza- 
tion in the diffusive regime to the results of Yu and 
Flatted based on standard drift-diffusion theory, we find 
that the field-dependent spin diffusion length L € s given 
by Eq. I|4.38|l agrees with the "up-stream" spin diffusion 
length L u given by Eq. (2.23b) of Ref. E provided the 
intrinsic spin diffusion length L of that reference is iden- 
tified with the spin diffusion length L s — ^/TT S of the 
present work. Then, by expressing the conductivity of 
the semiconductor in Eq. (3.5) of Ref. 19 (with the inter- 
face resistances set equal to zero) as a s = 2/3e 2 v e nil, we 
recognize the equivalence of that equation with our Eq. 
(|4.33|) in the diffusive regime. This, in turn, implies that 
the injected current and density spin polarizations of ei- 
ther work are formally identical. Numerical calculations 
have confirmed this result. 



3. Zero-bias limit 

We now consider the injected current spin polarization 
in the zero-bias limit, in which |ai| <C 1 and J(x) = J, 
i.e., re = 1. Here, the spin transport function A{x) is 
determined by Eq. i.e., A{x) oc exp(— x/L), so 

that Tj = 27/(1 + 7) = 7J and T n = 2/(1 + 7) = 7 „. 
From Eq. 14. 3311 . we then have 



Oil 



Pi 



X7j(l + Gi/ 7j g) 



(4.41) 



Combining this with Eq. 
current spin polarization 

Pj(xi) 



we find for the injected 



1 + GihjQ 



Pi 



(4.42) 



It is instructive to consider expression (|4.42|) in the dif- 
fusive and ballistic regimes. 



In the diffusive regime l/l s <C 1, we have 77 = 2y/l/l s 
and, therefore, 



Pj(xi) 



1 



1 + Gi/g 



Pi 



(4.43) 



Here, Qq = 2Q^Jljl s — <jq/L s , where Eq. (|3.7|l has been 
used, and (Jq = 2QI is the conductivity of the semiconduc- 
tor [see the remarks following Eq. (|2.71l) ]. The quantity 
Go is seen to be the semiconductor analogue of the fer- 
romagnet parameter G\ defined by Eq. I|4.4|) . Choosing 
Pi = 0.8 and adopting the values^ <j\ = 10 3 il _1 cm _1 , 
L s = 60 nm, <jq = 10 fJ^cm -1 , and L s — 2 /im, we 
have Gi/Go = 1.2 x 10 3 and hence Pj(xi) « 0.6 x 10" 3 . 
The large value of the ratio G\/Qq reflects the "conduc- 
tance mismatch" which appears to be the determining 
parameter of the injected spin polarization in the diffu- 
sive regime 

On the other hand, in the ballistic limit l/l s — > 00, we 
have 7,/ = 7 = 1, so that 



Pj(xi) 



1 



i + Gi/g 



Pi 



(4.44) 



Here, the Sharvin interface conductance Q takes the place 
of the quantity Qo in Eq. I|4.43|) . Assuming m* — 0.067m e 
and T = 300 K, we have Q = 3.2 x 10 11 frW 2 for 
n = 5 x 10 17 cm" 3 and Q = 0.64 x 10 10 Or 1 ^ 2 for 
n = 10 16 cm -3 . This results in Pj(xi) « 0.3 and 
0.8 x 10~ 2 , respectively. The first example, where the 
large Sharvin interface conductance entails a large in- 
jected spin polarization, is fictitious since the high dop- 
ing concentration needed to obtain an electron density of 
5 x 10 17 cm -3 (for example, in GaAs) would imply such 
small values of the momentum relaxation length / that 
ballistic transport is all but ruled out. Only semicon- 
ducting materials with unusually large mobilities would 
make this a realistic case. The second example with the 
lower electron density of 10 16 cm -3 would be more fa- 
vorable to a ballistic transport mechanism, but leads to 
a very small injected spin polarization; this confirms the 
conclusion of Kravchenko and Rashba 2 - stating that spin 
injection is suppressed even in the ballistic regime unless 
spin-selective interface resistances are introduced. 



V. CONCLUDING REMARKS 

We have developed a unified semiclassical theory 
of spin-polarized electron transport in heterostructures 
formed of a nondegenerate semiconductor and two ferro- 
magnetic contacts. In this theory, the spin polarization 
inside the semiconductor is obtained for a general trans- 
port mechanism that covers the whole range between the 
purely diffusive and purely ballistic mechanisms and is 
controlled by the momentum relaxation length of the 
electrons. 

The basis of the present work is provided by our previ- 
ously developed unified model of (spinless) electron trans- 
port in semiconductors, in which diffusive and ballistic 
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transport are combined in the concept of the thermobal- 
listic electron current. As a prerequisite to the exten- 
sion of the spinless unified model to spin-polarized trans- 
port, we have modified and completed its formulation in 
such a way that an unambiguous description of electron 
transport in terms of a uniquely defined thcrmoballistic 
chemical potential is achieved. From the chemical poten- 
tial, the unique thermoballistic current and density are 
obtained; numerical calculations show that, for typical 
parameter values, the thcrmoballistic current is close to 
the physical current. 

In order to treat spin-polarized transport in semicon- 
ductors within the unified description, we have intro- 
duced a thermoballistic spin-polarized current and a ther- 
moballistic spin-polarized density by allowing spin re- 
laxation to take place during the ballistic electron mo- 
tion. These are expressed in terms of a spin trans- 
port function which comprises in a compact form the 
information contained in the spin-resolved thermobal- 
listic chemical potentials. Using the balance equation 
that connects the thermoballistic spin-polarized current 
and density, we have derived an integral equation for 
the spin transport function, from which the latter can 
be calculated in terms of its values at the interfaces of 
the semiconductor with the contacts. The spin trans- 
port function determines, in conjunction with the spin- 
independent thermoballistic chemical potential, all spin- 
dependent quantities in the semiconductor, in particular, 
the position dependence of the current and density spin 
polarization. The spin polarization all across a ferro- 
magnet/semiconductor heterostructure is determined by 
making use of the continuity of the current spin polariza- 
tion at the contact-semiconductor interfaces and connect- 
ing the spin-resolved chemical potentials there. Thereby, 
a unified description of spin-polarized transport emerges 
that provides a basis for the systematic study of the in- 
terplay of spin relaxation and transport mechanism in 
heterostructures relevant to spintronic applications. 

To interpret the formalism developed here and to re- 
late it to previous, less general formulations, we have 
considered spin-polarized electron transport in a homo- 
geneous semiconductor without space charge, driven by 
an external electric field. Within a judicious approxima- 
tion, the integral equation for the spin transport function 
can then be reduced to a second-order differential equa- 
tion which generalizes the standard spin drift-diffusion 
equation to the case of arbitrary values of the ratio of 
momentum to spin relaxation length. In the zero-bias 
limit, the position dependence of the spin polarizations 
across a heterostructure is obtained in closed form. 

The generalized spin drift-diffusion equation 
has been used in calculations of the current 
spin polarization across a symmetric ferromag- 
net/semiconductor/ferromagnet heterostructure with 
material parameters in the range of interest for spin- 
tronic devices. The dependence on the transport 
mechanism in the semiconductor has been exhibited by 
varying the momentum relaxation length over several 



orders of magnitude. It was found that the ballistic 
regime favors sizeable (large) spin polarizations. The 
same picture emerges from calculations of the injected 
current spin polarization as a function of an applied 
electric field. While the field also serves to raise the 
polarization in the semiconductor, the main effect still 
is due to the variation of the momentum relaxation 
length, i.e., to the influence of the ballistic component 
of the transport mechanism. In order to exploit the 
potentiality of varying the transport mechanism with 
the aim to improve the efficiency of spintronic devices, 
the identification and design of novel semiconducting 
materials is called for. 

In the present work, emphasis has been placed on a 
careful elaboration of the formalism underlying the uni- 
fied description of spin-polarized electron transport in 
ferromagnet/semiconductor heterostructures. In the il- 
lustrative calculations, we have restricted ourselves to 
the simplest cases. In future work, applications of the 
present formalism will have to be based on the solution 
of the fundamental integral equation for the spin trans- 
port function in its general form. These should include 
the treatment of magnetic semiconducting samples (char- 
acterized, in the unified description, by nonzero values of 
the relaxed spin fraction excess) and of interface barriers, 
like Schottky or tunnel barriers (represented by appropri- 
ately chosen potential profiles). As to possible extensions 
of the theory, setting up a formalism for the treatment of 
degenerate semiconductors appears to have first priority. 

APPENDIX: UNIQUE THERMOBALLISTIC 
FUNCTIONS 

In this Appendix, we present details of the construc- 
tion of a unique thermoballistic chemical potential n(x), 
current J(x), and density n(x) in terms of the solutions 
Xi{x) and Xi{ x ) of Eqs. I|2.28|l and 1)2. 31|) . respectively. 

Evaluating expression (|2.12|) with the function J\(x) 
following from the solution of Eq. (|2.28|l [case (i)], with 
J 2 set equal to J\(x 2 ~)^ we obtain the thermoballistic 
current 

Ji(x) = wi(x 1 ,x 2 ;l) [Ji - Ji(x 2 )} 
f x dx' 

+ / - j -w 2 (x',x 2 ;l)[Ji(x')-J 1 (x 2 )} 
f X2 dx' 

+ / — w^xux'-^Ux-Jxix')] 

(A.l) 

and, similarly, for case (ii), the current 
J 2 {x) = w-i(xi,x 2 ;l) [J 2 {xi) - J 2 ] 

r x dx' 

+ / -j-w 2 (x',x 2 ;l)[J 2 (x')-J 2 } 
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r 2 dx' 

+ J —w 2 (x 1 ,x';l)[J 2 (x 1 )-J 2 (x')] 

+ [ T" /*' T" W3(x '' x " ; [j2(x,) " j2(2; " )] ' 

(A.2) 

which are not, in general, equal. This ambiguity is re- 
moved by introducing a unique thermoballistic current 
j( u >{x) as a superposition of the currents J\{x) and 
■Ma:), 

jM(x)=a 1 ^-J 1 (x)+a 2 ^-Mx) , (A.3) 

Jl ^2 

where 



■M 3 ^ ) 
J 2 



wi(xx,x 2 ;l) X2(xi 
a 

7 



2:2 dx' 

— w 2 (x',x 2 ;0x2(a:') • (A.ll) 



Employing Eqs. lfO|l . fO|l . and (|IT8)l - llA~TT|l . we ob- 
tain the coefhcients &i ; 2 as given by Eqs. I|2.35[) - H2.38[) . 
Now, introducing, in analogy to Eq. (|A.3|> . 



for X\ < x < x 2 , and, in addition, 



J {u \x x ) = a 1 ^-3 1 +a 2 ^-J 2 (x 1 ) 



■J 



J 



Ji 



(A.12) 



(A.13) 



1 f X2 

Ji,2 = / dxJ-xAx). (A.4) 

x 2 - xi J X1 

The current j( u \x) has to satisfy Eq. (|2~2U|) . 

\ r x 2 



dxJ iu) (x) = J 



(A.5) 



x 2 - xx J X1 

Similarly, Eq. (|2.24l) is to be replaced with 

jW(xf) = J {u \x 2 ) = kJ. (A.6) 

The conditions (|A.5|I and (| A. 6|> determine the coefficients 
ax and a 2 . We find from Eq. (|A.5(1 . using Eqs. (|A.3|) and 



ax + a 2 = 1 



(A.7) 



In order to apply condition l|A.6|) . we evaluate the cur- 
rents IjA.f I) and l|A.2|) at the ends of the sample, using 
Eqs. (|2~27jl and jZHjl , 



Jijxj) 
Ji 



= wx(xx,x 2 ;l) xi(x 2 ) 
f X2 dx' 

+ / —w 2 (xx,x';l)xi(x'), (A. 



J 2 (x+) 



= wi(xi,x 2 ;0 X2(a:i) 
f X2 dx' 

+ J —w 2 (xx,x;l)\x 2 (xx)-X2(x)} 



(A.9) 



) 

Jx 



wx(xx,x 2 ;l) xi{x 2 ) 



rX2 dx' 

+ I —w 2 (x',x 2 ;l)[xx(x 2 )~xi(x')] 



(A.10) 



J 



J w (x 2 ) =ax — Jx(x 2 ) + a 2 — J 2 
Jx J 2 



(A.14) 



we may write the unique thermoballistic current (jA.3Jl . 
using Eqs. I|A.1|) . I|A.2|) and (|A.12|) - (|A.14jl . in a symbolic 
form analogous to expression l|2.f 511 . 



J (M) (z) 



x dx[ 



dx 1 ' 



J^{x')-J^(x") 



W(x',x";l) 



(A.15) 



here, the values J^ u \xx, 2 ) are to be identified with their 
physical values ^71,2 in the contacts, 



J [u) {x ia ) =Jx, 2 • 



(A.I6) 



In line with the definition l|2.6(l of the current J(x) in 
terms of the chemical potential /i(x), we now define a 
unique chemical potential ^ u \x) via relation IjA.f 2|) by 



v„N, 



' e^J^\x) 



for xx < x < x 2 , where now 

V [u) {x ia ) = ^ 



(A.17) 



(A.18) 



The chemical potential fv u > (x) is the key quantity in the 
extended unified description of electron transport inside 
the sample. In terms of /i^ (x), the unique ballistic cur- 
rent across the interval [x',x"] appearing in the expres- 
sion for the thermoballistic current (|A.I5|) is given by 
Eq. H3J. 

For the explicit calculation of the chemical potential 
(i M (x), we use Eq. (fQ6|) in Eqs. l|A~T3)l and (|A~I4)l . 
We then find, with the help of Eqs. l(2~27jl and lj2~3U|) . 



1- ax^-\ Jx - a 2 ^-J 2 = a 2 Jx2(xx) , (A.19) 



Si j- Jx + (l-a 2 ^- ) J, = -a i J \ i ( .co ) 



(A.20) 
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By subtracting Eq. (|A.20(1 from Eq. HA.19(1 . we obtain 
Eq. I|2.41|) . On the other hand, adding Eqs. (|A.19f) and 
1|A.20|> results in 

ai- r Ji + a 2 - r J 2 = M+J2) 

,J\ J 2 2 

+ ^-[01X1(^2) -a 2 X2(xi)] ■ 

(A.21) 

Then, expressing the current J^ u \x) in terms of the 
quantities J7i,2 and the functions Xi,2(^) by combining 
Eq. ljA~T2)l with Eqs. l(2~27jl and (j!Hf])) . 

J {u) {x) = a^U.-.hxiix)} 

+ a 2 ^-[J 2 + J 2X2 (x)} 
J 2 

= 817 Ji+a 2 — J 2 
•h J 2 

- J[&ixi( x ) -02X2(2;)] , 



(A.22) 

we obtain Eq. (|2.39|l with X-{ x ) given by Eq. I|2.40|l . 
Using Eqs. I|2.39|l and 12.41(1 to eliminate the total current 
J as well as Eqs. (|A.17() to go over to the unique chemical 
potential p^ u '{x), we find 



"1 X-{x)~ 




'1 X-{x)~ 


[2 X „ 




[2 X 



(A.23) 

where x is defined by Eq. (|2.34|) . The corresponding 
thermoballistic current J^ u \x) and density rS u \x) are 
obtained by substituting expression (|A.23|I in Eq. (|2.1rj|) 
and H2.17fl . respectively. 

In the main body of the paper, we always deal with the 
unique chemical potential, current, and density, and omit 
the superscript u; we have already adhered to this con- 
vention when referring from the Appendix to the equa- 
tions of Sec. II. B. 
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